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Abstract 



Fungicide mixtures produced by the agrochemical industry often contain low risk fungi- 
cides, to which fungal pathogens are fully sensitive together with high risk fungicides 
known to be prone to fungicide resistance. Can these mixtures provide adequate dis- 
ease control while minimizing the risk for the development of resistance? We present 
a population dynamics model to address this question. We found that the fitness cost 
of resistance is a crucial parameter to determine the outcome of competition between 
the sensitive and resistant pathogen strains and to assess the usefulness of a mixture. 
If fitness costs are absent, then the use of the high risk fungicide in a mixture selects 
for resistance and the fungicide eventually becomes nonfunctional. If there is a cost of 
resistance, then an optimal ratio of fungicides in the mixture can be found, at which 
selection for resistance is expected to vanish and the level of disease control optimized. 



1. Introduction 



Fungal pathogens of plants pose a serious problem for sustainable food production [TJ. 
In order to protect crops from fungal diseases, farmers routinely apply fungicide sprays. 
However, in many cases plant pathogenic fungi develop resistance to one or more fungi- 
cides or fungicide classes, which makes the treatment ineffective and may lead to severe 
crop losses and lower profit for farmers. 

Fungicide resistance is a prime example of adaptation of a population to an environ- 
mental change, also known as evolutionary rescue j2j [3]. While global climate change 
is expected to result in a loss of biodiversity in natural ecosystems, evolutionary rescue 
is seen as a mechanism that may mitigate this loss. In the context of crop protection 
the point of view is quite the opposite: adaptation of pathogens of crop plants to the 
use of chemical disease control measures is seen as a major threat to food production. 
Therefore, any methods that diminish the adaptive ability of pathogens and decrease 
their diversity would be helpful to develop better crop protection strategies. This re- 
quires a detailed quantitative understanding of the dynamics of infection and the factors 
driving the emergence and development of fungicide resistance [I]. Despite the global 
importance and urgency of fungicide resistance, this problem has received relatively little 
theoretical consideration (see [51 El CH El E] and [I] for a comprehensive review). 

In recent years, agrochemical companies have begun marketing mixtures that contain 
fungicides with a low risk of resistance with fungicides that have a high risk of resistance. 
In extreme cases the high risk fungicide is no longer effective against some common 
pathogens because resistance has become widespread. For example, a large proportion of 
the European population of the important wheat pathogen Mycosphaerella graminicola 
(recently renamed Zymoseptoria tritici) [TOj [11] is resistant to strobilurin fungicides [12] 
but farmers continue to use the strobilurins in combination with other fungicides because 
strobilurins are thought to contribute additional physiological benefits or to protect the 
crop from other pathogens. 

Can the application of mixtures of low risk and high risk fungicides bring any advan- 
tage in terms of crop protection while minimizing further selection for resistance? Under 
which conditions does this type of mixture add anything to the treatment benefit, com- 
pared to when a low risk fungicide is applied alone? We addressed these questions using 
a simple population dynamics model of host-pathogen interaction based on a system of 
ordinary differential equations. We found that the fitness cost associated with resistance 
mutations is a crucial parameter, which governs the outcome of the competition between 
the sensitive and resistant pathogen strains. 

In the absence of a fitness cost, application of a high risk fungicide in a mixture 
inevitably and irreversibly selects for resistance, which makes its use in the mixture 
ineffective. On the other hand, if the resistance mechanism does confer a significant 
fitness cost to the pathogen, then application of the fungicide mixture may provide a 
stable advantage. In this case, an optimal ratio of the high risk fungicide in the mixture 
can be found below which the application of a high risk fungicide is beneficial, but 
at the same time does not select for resistance. The optimal ratio is determined by 
the cost of resistance, the degree of pharmacological interaction between the fungicides, 
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and depends on the degree of resistance for the case of partial resistance. Therefore, 
knowledge of the cost of resistance is necessary to determine whether the application of 
a given fungicide mixture will select for the resistant strain. 

A single point mutation associated with fungicide resistance often makes the pathogen 
completely insensitive to a fungicide, as is the case for the G143A mutation giving 
resistance to strobilurin fungicides in many fungal pathogens [131 HI]- In many other 
cases the resistance is partial and has properties of a quantitative trait, for example, 
resistance of Z. tritici and other fungi to azole fungicides [T5]. It was found that the 
resistance in Z. tritici results from accumulation of many mutations of the CYP51 gene, 
which encodes for the sterol 14a-demethylase, an enzyme involved in sterol biosynthesis 
and the target protein of azoles [T6l \T7\ ITS] . In order to account for these features, we 
have considered varying degrees of resistance and fitness costs in our model. 

Most of the previous modeling studies on the effect of fungicide mixtures on the 
development of the resistance were based on "density-independent" models, in which 
the pathogen population is assumed to grow exponentially over time and the effect of 
a finite host population is neglected (for example, [19], EDI EI])- However, the actual 
disease progress curves can usually be approximated by an exponential function only 
at the very beginning of an epidemic [22]. In a recent study [5j, a more realistic model 
was considered, in which the finiteness of the host tissue susceptible to infection was 
taken into account. However, in contrast to our study, resistance was assumed to bear 
no fitness costs for the pathogen. It was found that in the absence of fitness costs the 
use of fungicide mixtures delays the development of resistance [5] . This conclusion is in 
agreement with our results. In the absence of fitness costs the application of the high 
risk fungicide gives a selective advantage to the resistant pathogen strain. But the rate 
at which the frequency of the resistant strain increases is proportional to the difference 
in fitness between the resistant and sensitive pathogen strains: the larger the difference, 
the faster is the selection. The sensitive pathogen strain is suppressed in the same way 
by the high risk and the low risk fungicide, therefore, its fitness is the same when the 
high risk fungicide is applied alone and when a mixture of the two fungicides is applied. 
In contrast, the fitness of the resistant strain is diminished to a greater extent by the 
mixture than by the high risk fungicide alone. Thus, the difference in fitness between 
the two strains is smaller and the selection for resistance is delayed when the mixture is 
applied (see Appendix lA.4l in the Electronic Supplementary Material for more details). 
Here we focus on finding conditions under which the selection for the resistant pathogen 
strain is prevented rather than delayed by using fungicide mixtures. 

The problem of combining chemical biocides in order to delay or prevent the develop- 
ment of resistance also appears in other contexts, such as resistance of agricultural weeds 
to herbicides [23], insect pests to insecticides [21], and human pathogens and diseases 
including HIV [22], bacterial infections [2E], malaria [2Z] and cancer [2S]- Accordingly, 
lessons from fungicide combinations may well apply to the problem of biocide resistance 
in many other contexts. 
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2. Model and assumptions 



In order to investigate the effect of fungicide mixtures on selection for fungicide resis- 
tance, we use a deterministic mathematical model of susceptible-infected dynamics (see 
Fig.HD 

^L = r(K-H-I s - J r ) - 6 [I s + (1 - p T )I T ] H, (1) 
at 

^ = bHI s -[l + e s (C,r B )] f iI s , (2) 

^ = 6(1 - Pr )HI r - [1 + S r (C, r B )] filr- (3) 

The model has three compartments: susceptible hosts H, hosts infected by a sensitive 
pathogen strain I s , hosts infected by a resistant pathogen strain 7 r ; and is similar to 
the models described by [HE]. The subscript "s" stands for the sensitive strain and the 
subscript "r" stands for the resistant strain. The quantities H, I s and I T , represent the 
total amount of the corresponding host tissue within one field, which could be leaves, 
stems or grain tissue, depending on the specific host-pathogen interaction. Susceptible 
hosts H grow with the rate r. Their growth is limited by the "carrying capacity" K, 
which may imply limited space or nutrients. Furthermore, susceptible hosts may be 
infected by the sensitive pathogen strain and transformed into infected hosts in the 
compartment I s with the transmission rate b. Susceptible hosts may also be infected by 
the resistant pathogen strain and transformed into infected hosts in the compartment 
I T . In this case, the resistant strain suffers a fitness cost p T which affects its transmission 
rate such that it becomes equal to 6(1 — p r ). The corresponding terms in Eqs. ([I])-© 
are proportional to the amount of the available susceptible tissue H and to the amount 
of the infected tissue I s or I T . Infected host tissue loses its infectivity at a rate [i due to 
pathogen death. 

Since our description is deterministic we do not take into account the emergence of 
new resistance mutations but assume that the resistant pathogen strain (J r ) is already 
present in the population. Therefore, when "selection for resistance" is discussed below, 
we refer to the process of winning the competition by this given resistant strain due 
to its higher fitness with respect to the sensitive strain in the presence of fungicide 
treatment. Emergence of new resistance mutations is a different problem, which goes 
beyond the scope of our study and requires stochastic simulation methods. We do not 
consider the possibility of double resistance in the model, but by preventing selection 
for single resistance as described here, one would also diminish the probability of the 
emergence of double resistance for both sexually and asexually reproducing pathogens 
(see Appendix lA.7p . 

We consider two fungicides A and B. The fungicide A is the high risk fungicide, to 
which the resistant pathogen strain exhibits a variable degree of resistance. However, 
the sensitive strain is fully sensitive to fungicide A. The fungicide B is the low risk 
fungicide, i. e. both pathogen strains are fully sensitive to it. We compare the effect of 
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the fungicide A applied alone, fungicide B applied alone and the effect of their mixture 
with an arbitrary proportion. We assume that the fungicides will increase the pathogen 
death rate [see the expression in square brackets in Eq. §2§, Eq. (J3J)]. Whether only one 
or a mixture of two fungicides is applied, the general form of the fungicide action is 
described by 

e(C) = K-^- (4) 

^ + ^50 

This function grows with the fungicide concentration C and saturates at a value k^, the 
maximum killing rate. The parameter C50 is the fungicide concentration at which half 
of the maximum effect is achieved; it determines how fast the function e(C) grows. We 
assume that the action of the fungicide A when applied alone is described by Eq. 01]): 
e(Ca) = k-^xCx/{Cx + C 50 a) an d the action of the fungicide B when applied alone is 
also given by Eq. (j4j): s(Cb) = k^C-s/i^C-s + Csob)- For simplicity we assume that 
&kA = &kB = &k and C50A = C50B = Cso- If the parameters C50A and C50B are different, 
the concentration axis can be rescaled so that they become the same. A more general 
case, where the maximum killing rates and k^B are different can also be considered, 
but leads to slightly more complicated expressions. The function e(C) in Eq. 01]) is a 
special case of the sigmoid Emax model [29] with the Hill coefficient equal to one, which 
is often used in the antibiotic [30J and fungicide resistance literature [6]. 

When a mixture of two fungicides is applied, it increases the death rate of the sensitive 
pathogen according to 

C 

e s {C,r B ) = h , (5) 

and the death rate of the resistant pathogen is increased by 

C 

£ r (C, r B ) = A: k . (6) 

^ + C 50 /7r 

Here C = C*a + Cb, where C*a is the concentration of the fungicide A and Cb is the 
concentration of the fungicide B, r B = Cb/C is the proportion of the fungicide B in the 
mixture and 

7s = 1 + uy/r B (l-r B ), (7) 



7 r = a(l - r B ) + r B + u\/ar B (l - r B ) (8) 

are the parameters which modify C50 due to pharmacological interaction and partial re- 
sistance. The Eqs. ([5]), ([6]) are obtained from more general Eqs. (1A.4I) . (1A.5j) by changing 
the variables from Ca, C b to C, r B and assuming that the sensitive pathogen strain is 
fully sensitive to both fungicides and the resistant strain can have a varying degree of 
sensitivity a to the fungicide A, but is fully sensitive to the fungicide B. The degree of 
sensitivity a enters as the coefficient before (1 — r B ) in Eq. (jSJ) and lies between zero and 
one. At a = the pathogen is fully resistant to fungicide A, and at a = 1 the pathogen 
is fully sensitive to it. 
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The degree of pharmacological interaction is characterized by the parameter u. At 
u = the fungicides do not interact, u > represents synergy and u < is the case of 
antagonism. (See Appendix lA. 1 1 for the discussion of pharmacological interaction.) Note 
that the term responsible for the fungicide interaction in Eqs. (EJ), (JH]) is proportional to 
^/tb(1 — tb), which is zero when only one of the fungicides is applied, i. e. at re = or 
r-Q = 1, and reaches a maximum when the two fungicides are mixed in equal proportions, 
i.e. r B = 1/2. 

We consider one growing season and assume the host-pathogen equilibrium to be 
reached during the season. That corresponds to the time-dependent solution of Eqs. ([T])- 
fl3]) reaching its stable steady state. When saturation of the empirical disease progress 
curve occurs at a value below 100 %, this is a good indication that the equilibrium is 
reached during one season. Whether this is the case depends on both the host-pathogen 
combination and the environmental conditions. While in some cases the disease progress 
curve saturates (for example, [HUES]), in other cases it does not (for example, [3"2"j I2"2"]). 

We assume that the fungicide concentration is constant over time and neglect the 
spatial dependencies of the variables H, I s and I r and all other parameters. The latent 
phase of infection, which could be considerable for some pathogens (for example, Z. tritici 
has a latent period of about 12 days), is also neglected. Since we neglect mutation, 
migration and spatial heterogeneity, the two pathogen strains cannot co-exist in the 
long-term. 

Such competitive exclusion is observed in the field in some cases, but in other cases 
resistant and sensitive pathogen strains are found to coexist. For example, a gradual 
increase in the frequency of fungicide-resistant pathogen strains up to almost 100 % 
was observed in field experiments in plots treated with fungicides, but a stable co- 
existence of resistant and pathogen strains developed in untreated plots [33] • In a field 
monitoring study [M] a fungicide-resistant pathogen strain emerged under fungicide 
treatment and gradually dominated the pathogen population over the course of five years. 
A modeling study [7] has shown that spatial inhomogeneity in the form of imperfect 
fungicide coverage makes the stable co-existence of resistant and sensitive pathogen 
strains possible if they have similar values of the basic reproductive number. 

The basic reproductive number (denoted as Ro) is often used in epidemiology as a 
measure of transmission fitness of infectious pathogens [35]. It is defined as the expected 
number of secondary infections resulting from a single infected individual introduced into 
a susceptible population. At R > 1 the infection can spread over the population, while 
at Rq < 1 it cannot and the epidemic dies out. 

The equilibrium stability analysis of the model Eqs. (JTJ)- (|3j) shows that the relationship 
between the basic reproductive number of the sensitive strain 



and the basic reproductive number of the resistant strain 



-Ror — 77— — TrT~\T~ (10) 
[1 + e r (C, r B )]/i 

determines the long-term outcome of the epidemic. The sensitive strain wins the com- 
petition and dominates the pathogen population if Rq s > 1, such that it can survive in 
the absence of the resistant strain, and Rq s > Rq t , such that is has a selective advantage 
over the resistant strain. The second inequality is equivalent to 

< (Pr/h+ „ , g , )/(l-Pr). (11) 



C + C 50 / 7 s V C + C 50 /7 r 

Similarly, the resistant strain wins the competition and dominates the population if 
Rq v > 1 and i?or - > Rqs- 

We assumed here that resistance cost decreases transmission rate b of the resistant 
pathogen strain and the fungicide increases the pathogen death rate \i. We performed the 
same equilibrium analysis when the effect of the resistance cost and the fungicide enter 
the model in other ways and obtained qualitatively the same results (see Appendix lA.5p 

The simplicity of the model allows us to obtain all the results analytically. We de- 
termined explicit mathematical relationships between the quantities of interest, which 
enabled us to study the effects over the whole range of parameters. 



3. Results 

Understanding the factors determining the outcome of the competition between sensitive 
and resistant pathogen strains is important for developing an optimal treatment regime. 
Hence, we first investigate the parameter ranges over which resistant or sensitive strains 
dominate the pathogen population for the case of fungicides A and B applied individually 
and for the mixture of the two fungicides (Sec. l3.lj) . Then, we consider the optimal 
proportion of fungicides to include in a mixture in Sec. 13.21 and the benefit of fungicide 
treatment in Sec. 13. 31 Finally, we take into account possible pharmacological interactions 
between fungicides and consider the effect of partial resistance (Sec. l3.4[ 13. 5p . 



3.1. Selection for resistance 

The ranges of fungicide concentration and cost of resistance at which the sensitive (white) 
or resistant (dark red) pathogen strain is favored by selection are shown in Fig. [2] If a low 
risk fungicide is applied alone, the sensitive strain has a selective advantage across the 
whole parameter range in Fig.^a). When only a high risk fungicide is applied [Fig.^b)] 
the resistance dominates if the fitness cost is low p r < k^/{k^ + 1) and at a fungicide 
concentration higher than a threshold value which increases with the fitness cost [see 
Eq. (1A.22P ], If the fitness cost exceeds k\/{k^-\- 1) [dotted vertical line in Fig.E^b)], then 
the sensitive strain dominates at any fungicide concentration. 
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The Fig.[2]^c) shows the outcome when the two fungicides are mixed at equal concen- 
trations: here the picture changes considerably. Now the borderline between the area 
where the sensitive strain dominates at any fungicide concentration and the area where 
resistance may dominate, is shifted to smaller values. Moreover, to the left of the bor- 
derline, an area where the sensitive strain dominates appears at concentrations above a 
threshold value. 

Importantly, without a fitness cost (p r = 0) the resistant strain becomes favored by 
selection and will eventually dominate the population whenever the high risk fungicide 
is applied, alone or in combination with the low risk fungicide [Fig.|2^b,c)]. 

Interestingly, the value of the borderline [vertical dotted line in Fig.^c)] depends on 
the ratio of the fungicides in the mixture. This allows one to determine an optimal 
proportion of fungicides to use in the mixture (see Sec. l3.2|) . 



3.2. Optimal proportion of fungicides in a mixture 

It is highly desirable to keep existing fungicides effective for as long as possible. There- 
fore, an optimal mixture contains the largest proportion of the high risk fungicide, at 
which (i) the resistant pathogen strain is not selected and (ii) an adequate level of disease 
control is achieved. In order to fulfill both of these objectives, it is necessary for the 
fitness cost of resistance to be larger than a threshold value p r > p r b [see Eq. flA.13j) in 



Appendix]. This means that selection is operating to the right of the vertical dashed line 
in Fig.|2]^c), where the sensitive strain is favored and the resistant strain is eliminated 
at any fungicide concentration. 

The threshold value p r b depends on the proportion of fungicides in the mixture. 
Adding more of the low risk fungicide, while keeping the total concentration C con- 
stant, shifts the threshold to the left. This diminishes the range of the values for fitness 
cost at which the resistant strain is competitively superior. On the other hand, adding 
less of the low risk fungicide, while again keeping C constant, shifts the border to the 
right, increasing the parameter range over which the resistant strain is favored. 

Therefore, at a given fitness cost p r , one can adjust the fungicide ratio r-Q in such a 
way that p T > p r b- This is shown in Fig.|3j the curve shows the critical proportion of the 
low risk fungicide tbc, above which the treatment with a mixture is effective in the long 
run, meaning no selection for resistance occurs at any total fungicide concentration C. 
One can see from Fig. [3] that if the resistance cost is zero (p r = 0), then for the treatment 
to be effective the fungicide A should not be added at all, i.e. r-g = 1. As the fitness 
cost grows, the value of tb c decreases, giving the possibility to use a larger proportion 
of the high risk fungicide without selecting for resistance. 

So far we have shown how choosing an optimal proportion of fungicides in the mixture 
prevents selection for resistance. Now, we will consider how to achieve an adequate level 
of disease control. In order to do this, we need to quantify the benefit of fungicide 
treatment. 
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3.3. Treatment benefit 



The yield of cereal crops is usually assumed to be proportional to the healthy green 
leaf area, which corresponds in our model to the amount of susceptible hosts H(t). 
Accordingly, we quantify the benefit of the fungicide treatment as the ratio between 
the amount of susceptible hosts H(t) when both the disease and treatment are present 
and its value H nd (t) in the absence of disease: B(t) = H{t)/H nA {t). Hence, B{t) = 1 
corresponds to a perfect treatment, which eradicates the disease completely and the 
treatment benefit of zero corresponds to a situation where all susceptible hosts are 
infected by disease. This quantity characterizes the effectiveness of the chemical disease 
control. 

The treatment benefit at equilibrium as a function of the fitness cost and the fungi- 
cide concentration is shown in Fig.|H (The corresponding equations are given in Ap- 
pendix lA~3l ) The range of fitness costs and fungicide concentrations is the same as in 
Fig. El The upper graph shows the outcome when a low risk fungicide is applied alone. 
Here, the treatment benefit does not depend on the cost of resistance, but increases 
monotonously with the fungicide concentration, since the sensitive strain is favored by 
selection here. 

The dependency changes when a high risk fungicide is applied alone (middle graph): 
the benefit increases with the fungicide concentration and does not depend on the cost 
of resistance only in the parameter range where the sensitive strain is favored (to the 
right from the black solid curve). In the range where resistance dominates, the treat- 
ment benefit depends only on the cost of resistance and is independent of the fungicide 
concentration, as we assume resistant mutants are fully protected from the fungicide A. 
Hence, increasing its concentration does not affect the amount of disease. 

In the case of a mixture of a high risk and a low risk fungicide, the parameter range 
in which the treatment benefit depends only on the fungicide concentration and is inde- 
pendent of the cost of resistance becomes larger. In the area where resistance is favored 
(to the left from the black solid curve) the treatment benefit increases with both the 
cost of resistance and the fungicide concentration, since the fungicide B in the mixture 
suppresses the resistant strain. 

Both fungicides in the mixture contribute to the suppression of the disease only in 
the region of Fig.H^c), where the sensitive strain dominates (to the right from the black 
solid curve). Moreover, to achieve the desired treatment benefit by adjusting the total 
fungicide concentration C without the danger of selecting for resistance, the fitness cost 
should exceed its threshold value p r > p r b [to the right from the dashed vertical line 
in Fig.lHc)]. As we have shown above, the threshold value of the resistance cost p r b 
depends on the proportion of fungicides in a mixture. Thus by first choosing an optimal 
proportion of fungicides r B as suggested in Sec. l3.2l and then choosing the total fungicide 
concentration such that the treatment benefit reaches a high enough value, one can avoid 
selection for resistance and at the same time reach the desired level of disease control. 



10 



3.4. The effect of pharmacological interaction between fungicides 

Synergistic interactions between fungicides makes their combined effect greater than 
expected with purely additive interactions. The sensitive pathogen strain is suppressed 
more effectively by a synergistic mixture, while the resistant strain is not affected by 
the interaction (in case of full resistance a = 0). This increases the range of fitness 
costs over which resistance has a selective advantage [the dashed vertical line in Fig.[2^c) 
shifts to the right]. Consequently, the critical proportion of the low risk fungicide in 
the mixture r^ c , above which the resistant mutants are eliminated increases [dotted 
curve in Fig.E^a)]. In contrast, an antagonistic mixture suppresses the sensitive strain 
less effectively than either fungicide used alone. Hence the range of fitness costs over 
which resistance dominates becomes smaller and the ratio r-Q C decreases [dashed curve in 
Fig.[3(a)]. This result is in agreement with studies on drug interactions in the context of 
antibiotic resistance, where antagonistic drug combinations were found to select against 
resistant bacterial strains [36J. 

3.5. The effect of partial fungicide resistance 

Consider the situation when the resistant pathogen strain is not fully protected from 
the high risk fungicide, but exhibits a partial resistance (0 < a < 1). In this case, 
the fungicide mixture is more effective in suppressing the resistant strain than in the 
case of full resistance (a = 0) considered above. Therefore, one needs less of the low 
risk fungicide in the mixture to reach the conditions where resistance is eliminated 
by selection: the critical proportion of the low risk fungicide in the mixture decreases 
with the degree of sensitivity a in Fig.H^b). Also, in Fig.H^a) the dependency of the 
critical ratio of the fungicide B in the mixture for partial resistance (light grey curve) 
lies below the one at perfect resistance and reaches zero at a much smaller value. Thus, 
knowledge of the degree of resistance is crucial for determining an appropriate proportion 
of fungicides in the mixture. 

4. Discussion and conclusions 

The above results demonstrate that if fungicide resistance comes without a fitness cost, 
application of fungicides prone to resistance (high risk fungicides) in a mixture with 
fungicides still free from resistance (low risk fungicides) will select for resistance. Con- 
sequently, as soon as the equilibrium is reached the resistant strain will dominate the 
pathogen population and the sensitive strain will be eliminated. Because of this, the 
high risk fungicide will not affect the amount of disease and only the low risk fungicide 
component of the mixture will be acting against disease. Hence, the high risk fungicide 
becomes nonfunctional in the mixture and using the low risk fungicide alone would have 
the same effect at a lower cost. 

We extensively searched the available literature on fitness costs in different fungal 
pathogens of crop plants. A few studies inferred significant fitness costs from field 
monitoring (see for example, [37] and references in [38] )• But the findings in these studies 
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could also be due to other factors, including immigration of sensitive isolates, selection 
for other traits linked to resistance mutations or genetic drift [38]. Though relatively 
few carefully controlled experiments have been conducted, the majority indicate that 
fitness costs associated with fungicide resistance are either low (for example, [23110]) or 
absent (for example, [HJ H2])- Moreover, even if fitness costs are measured in resistant 
mutants generated in the lab [39], they might not represent those selected in the field, 
because compensatory mutations improving pathogen fitness are likely to occur under 
field conditions [35] • Thus the use of fungicide mixtures that include high-risk or already 
failed fungicides will select for resistance for most plant pathogens and cannot usually 
be justified according to the results from our model and current knowledge on fitness 
costs. But in the rare cases where fitness costs are significant (for example [12]), high 
risk fungicides in a mixture with low risk fungicides are expected to contribute to the 
suppression of disease and not to select for resistance if used at a low enough proportion. 

However, current knowledge of fitness costs is insufficient, since they may manifest in 
many different ways (for example, by impairing the ability of a pathogen to overseason 
[H]) and depend on environmental conditions. We hope this study will stimulate further 
experimental investigations to better characterize them and expect that significant costs 
will be found in some cases. In these cases, high risk fungicides can be used effectively 
for an extended period of time. An optimal proportion of the high risk fungicide in a 
mixture with the low risk fungicide can be determined that contains as much as possible 
of the high risk fungicide, but still does not select for resistance while providing adequate 
disease control (see Box[T]). If a mixture with the optimum proportion is applied, then 
the rise of the resistant strain is prevented for an unlimited time. Under the scenario 
where the resistant pathogen strain is fully protected from the high risk fungicide, the 
fungicides do not interact and the fitness cost of resistance decreases transmission rate 
b by 5 percent, then up to 20 percent of the high risk fungicide can used in the mixture 
without selecting for resistance (Fig. [3]). 

While it was previously discussed [15] that alternation of high risk and low risk fungi- 
cides might be a useful tactic for disease control in the presence of a fitness cost, we 
have shown that a mixture of these fungicides in an appropriate proportion can provide 
adequate disease control without selecting for resistance. Mixtures offer an advantage 
compared to alternation because there is no need to delay the application of the high 
risk fungicide and the resistant strain does not rise to high frequencies, which lowers the 
risk of its further spread (see Appendix lA.6p . 

This analysis applies only for a known resistant strain or mutation. If a new resistance 
mutation emerges with a different fitness cost, then it should also be taken into account 
when determining the optimal proportion of fungicides in the mixture. Alternatively, if a 
compensatory mutation emerges which strongly diminishes the fitness cost of resistance, 
then the resistant strain carrying the compensatory mutation may rise in frequency and 
eventually outcompete the sensitive strain. 

Our model considered only one pathogen that has a resistant strain, but there can 
be several different pathogens with resistant strains infecting a given host population. 
To minimize costs farmers prefer to use the same fungicide mixture to control as many 
pathogens as possible. Our approach could be generalized for this case: a proportion of 
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In order to avoid selection for resistance while providing adequate disease control, one 
should choose the fungicide ratio r-g and the total concentration C in the following way: 

1. measure the pharmacological properties of both fungicides to determine k^, and 
C50. 

2. determine the degree of pharmacological interaction u between fungicides A and 
B. 

3. measure the fitness cost of resistance p r 

4. choose the proportion of the fungicide B above the threshold: r-Q > rsc, such that 
the resistance is not favored by selection at any total fungicide concentration C. 

5. choose the total fungicide concentration, according to the desired level of disease 
control (see Fig.H]). 

Box 1: How to determine an optimal mixture of fungicides. 

fungicides in the mixture can be determined at which resistant strains of all pathogens 
are not favored by selection. 

The optimal proportion of fungicides in the mixture depends on the fitness cost of 
resistance, the degree of resistance, dose-response parameters of the fungicides and the 
degree of their pharmacological interaction. All these factors affect the optimal fungicide 
proportion which may vary from one host-pathogen combination to another. Hence, the 
end user should be able to mix fungicides in a controlled way (e. g. a tank mixture) 
in order to control a single pathogen or a group of pathogens common in the region. 
Thus, we agree with [3] that a tank mixture has an advantage compared to a formulated 
mixture (a mixture produced by the fungicide manufacturer). 

The theoretical results presented here need to be tested experimentally. We suggest 
to measure the amount of disease and the frequency of resistance as functions of time 
at different proportions of the high and low risk fungicides in the mixture. In this way, 
the existence of an optimal proportion of fungicides in the mixture could be tested em- 
pirically in cases where significant fitness costs are associated with fungicide resistance. 
The model prediction that the mixture will be ineffective if fitness costs are absent can 
also be tested in this way. 
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Figure 1: Scheme of the model in Eqs. ([I])-©. 
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Figure 2: Outcomes of the competition between the sensitive and resistant pathogen 
strains depending on the fitness cost of resistance p r and the fungicide con- 
centration C when treated with a single fungicide B at (Cb = C, panel (a)), 
a single fungicide A (Cb = C, panel (b)) and the combination of fungicides 
A and B (Ca = Cb = C/2, panel (c)). The range of the total fungicide con- 
centration C and the fitness cost of resistance p r , in which resistant strain is 
favored is shown in dark red. The range where selection favors the sensitive 
strain is shown in white. The vertical dashed line and the solid black curve in 
panel (b) are plotted according to Eq. (1A.21I) and Eq. flA.22j) in Appendix lA.2| 
respectively. The vertical dashed line and the solid black line in panel (c) are 
drawn according to Eq. (IA. 13[) and Eq. (1A.14I) in Appendix lA.il respectively, at 
7 S = 1, 7 r = 1/2. Fungicides are assumed to have zero interaction (u = 0) and 
the resistant strain is assumed to be fully protected from fungicide A (a = 0), 
the fungicide dose- response parameters are = 1, C50 = 1. 
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resistance cost p Y 

Figure 3: The critical proportion r-Q C of fungicide B (low risk fungicide) in the mixture, 
above which there is no selection for the resistant strain at any total fungicide 
concentration C, plotted according to Eq. (IA.17[) as a function of the resistance 
cost p r , assuming no pharmacological interaction (u = 0) and full resistance 
(a = 0). 
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Figure 4: Treatment benefit as a function of fungicide concentration C and fitness cost 
of resistance p T , plotted according to Eq. (1A.25j) in panel (a), according to 
Eq. (1A.26j) in panel (b) and according to Eq. f lA.27j) in panel (b). Treatment 
with fungicide B is shown in panel (a). Treatment with fungicide A is shown 
in panel (b). Treatment with a mixture of A and B at equal concentrations 
(r-Q = 1/2) is shown in panel (c). The solid and dashed black in panels (b) and 
(c) are the same as in Fig.[2j Fungicides are assumed to have zero interaction 
(u = 0) and the resistant strain is assumed to be fully protected from fungicide 
A (a = 0). The maximum killing rate — 1, the basic reproductive number 
of the sensitive strain without fungicide treatment Rq s (C = 0) = bKj fi = 2. 
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Figure 5: The effect of pharmacological interaction and partial resistance on tbc, the 
critical ratio of the fungicide B. r^ c is plotted as a function of the fitness cost of 
resistance p r (left panel), according to Eq. flA. 13[) for the case of no interaction 
between the fungicides u = (solid, the same as the curve in Fig.[3]), synergy 
u = 0.9 (dotted), and antagonism u = —0.9 (dashed) for the case of perfect 
resistance a — 1. The case of partial resistance at no interaction (a = 0.5, 
u = 0) is shown as a light grey curve. r-Q C is shown as a function of the degree 
of fungicide sensitivity a at p T = 0.05 (solid) and p T — 0.1 (dash-dotted) also 
according to Eq. (I A . 1 3 [) in the right panel. 
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A. Electronic Supplementary Material 



A.l. Model equations 

In order to elucidate the assumptions we made in Sec.El we consider a more general 
system of equations, which describes the change in time of the same quantities as in 
Eqs. ([I])-©: the amount of susceptible host tissue H, the amount of host tissue infected 
with the sensitive pathogen strain I s and the amount of host tissue infected with the 
resistant pathogen strain I T 

^L = r(K-H-I s -I T )-b [(1 - p s )I s + (1 - p r )I T ] H, (A.l) 

^ = 6(1 - Ps )HI s - [1 + e s (C A , C B )] iiI s , (A.2) 

^ = 6(1 - Pl )HI r - [1 + e T (C A , C B )] fil t , (A.3) 

where, the function e s (C A ,C B ) describes the effect of the application of the mixture 
fungicides A and B with concentrations C A and C B on the death rate of the sensitive 
pathogen strain and the function e v (C A ,C B ) describes the effect of this mixture on the 
death rate of the resistant strain: 

e s {C A ,C B ) = k 



£r(C A ,C B ) = k k 



a s , A C A + « s ,bCb + C 50 / [l + u^a StA a StB C A C B / (a s , A C A + a s>B C B ) 

(A.4) 

Qr.A^A + a T , B C B 

o^aCa + a r ,BC B + C 50 / [l + u^Ja^ A a^ B C A C B / (a r , A C A + a I>B C B )] 

(A.5) 



The parameters a s , A , a SjB , a T , A and a VjB characterize the degree of sensitivity of each of 
the two pathogen strains (index "s" for the sensitive strain, index "r" for the resistant 
strain) to each of the two fungicides A and B. Their values are between zero and one. 
In this general case both pathogen strains are partially resistant to both fungicides. 

The parameter C50 in Eqs. (1A.4[) . flA.5[) is modified due to pharmacological interaction 
between fungicides characterized by the degree of interaction u. At u = fungicides 
do not interact, u > represents synergy and u < corresponds to antagonism. [We 
restrict our consideration to u > —1, since otherwise the term in the square brackets 
of Eqs. (IA.4p . (1A.5P may become negative, which makes no sense.] This way to define 
pharmacological interaction between compounds is called "Loewe additivity" or "concen- 
tration addition" in the literature [451 HZ]. m this approach an interaction of a compound 
with itself is set by definition to be additive (zero interaction). For example, when the 
fungicide A is mixed with itself, the resulting sham mixture is neither synergistic, nor 
antagonistic but has zero interaction. An alternative way to define pharmacological 
interaction assumes that the two compounds have independent modes of action and is 
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called "Bliss independence" [48] or Abbott's formula [49]. However, in this definition 
a compound can have a pharmacological interaction with itself, i. e. be synergistic or 
antagonistic. 

There are several ways to introduce a deviation from the zero interaction regime, in 
which usually an interaction term is added to the isobologram equation [16]. We have 
chosen a specific form of the interaction term, which is proportional to the square root 
of the product of the concentrations of the two compounds [Eq. (28) in [46J]. This form 
allows for a simple analytical expression of the effect of the combination in Eqs. (1A.4[) . 

(ED- 

We assume that the cost of resistance decreases the transmission rate b by a fixed 
amount p s for the sensitive strain and by p T for the resistant strain in Eqs. f lA.lj) - flA.3j) . 



We restrict our consideration here to the case when the "sensitive" pathogen strain 
is fully sensitive to both fungicides (a Sj A = a s ,B = 1) and the "resistant" strain can 
have varying degrees of resistance to the fungicide A (a r) A = a, < a < 1), but is 
fully sensitive to the fungicide B (a T}B = 1). Therefore, the cost of resistance for the 
sensitive strain is zero p s = 0. Then, the fungicide dose-response functions become 
simpler [Eqs. (jHJ)]. 

In order to determine the range of fitness costs p T and fungicide concentrations C, 
over which the sensitive strain is favored by selection we perform the linear stability 
analysis of the fixed points of the system Eqs. (lA.lj) - (]A.3|) . which is equivalent to solving 
the inequalities Ro s > 1, -Ros > -Ror, where 

R ° s = [l + e s (C A ,C B )]^ (A ' 6) 
is the basic reproductive number of the sensitive strain and 

7? &(1 -ft)* ( , 7) 

is the basic reproductive number of the resistant strain. 

We consider then the inequality Ro s > Ro r , which is equivalent to 

< (Pr/h + „ , Z , ) /(l - Pr), (A.8) 



c + c 5 ohs V c + c 50 /i 

where 



7s 



l + uy/r B (l-r B ), (A.9) 



7 r = a(l - r B ) + r B + u v / ar B (l - r B ) (A. 10) 

and r B = C B /C is the proportion of the funcigide B in the mixture, C = C\ + C B . 
The inequality flA.8j) holds at 

p r < p rb , for (C < C b i or C > C h2 ) (A. 11) 
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or at 



Here, 



Px > Prb, for any value of C. (A. 12) 



C&1,2 

where 



&k(7s - 7r) (is + 7r + 7s^k - 2A/7r7s(l + &k)) 

( 7r -7s(l + £; k )r 



Prb = —, ; ; 2 -, (A. 13) 



27s7rPr(l + fck) L 



7s(p r -A; k + / 1 ; k p r )+7 r (/ 1 ; k + p r )±v A D , (A.14) 



D = 7s (^k(Pr " 1) + Pr) + 7r (*k + Pl f + 2 7s7r (^(p r - 1) - p\ - k k p 2 ) . (A.15) 

According to the inequality (IA.12I) . if the fitness cost of resistance is larger than a 
threshold value given by Eq. (IA.13I) . the sensitive strain has a selective advantage and 
the resistant strain is eliminated from the population at any fungicide concentration 
C > 0. 

For the case of no interaction between fungicides (u = 0) and perfect resistance (a = 0) 
the Eq. (1A. 13[) is simplified: 



(1 - r B )(l + k k + r B - 2y/r B (l + k k )) 

Prb = h — — . A.16 

(1 + k k - r B y 

We then solve the inequality p > p r b with respect to re and find that it is fulfilled at 
r B > r Bc, where 



P 2 T + k 2 k (l - p T ) + k kPr ( Pr - 2y/(l + k k )(l - Pr )) 

rBc = IhT^ (A - 17) 

It represents the critical proportion of the fungicide B in the mixture above which the 
resistant strain is not favored by selection (Fig.[3]). 



A. 2. Selection for resistance at no interaction between fungicides 

When only the high risk fungicide (fungicide A) is applied with the concentration C*a, 
we set r B = in Eq. (IA.90 and Eq. (lA.lOj) to obtain j s = 1, j T = a. Then, the following 
expressions are obtained for the threshold value of the resistance cost from Eq. flA.13j) 



k k (l - a) (l + a + kk - 2y/a(l + k k ) 

Prb = ±- -g '-, (A.18) 

{a -fa- 1) 
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and the fungicide concentration from Eq. flA.14j) 



c. 



50 



61,2 



2ap r (l + k k 



(p r - k k (l - p r )) + a(k k + p r ) ± VD 



where 



D = a 2 (k k + p r ) 2 - a (kl(l - p r l) + p 2 + Kpf) + p r - k k (l - p T ). 



(A.19) 



(A.20) 



In the simpler case of full resistance we take the limit a — > 0. Then, by taking this limit 
in Eq. (1A. 18[) . Eq. (1A. 19|) and Eq. (1A.20I) we obtain for the threshold values of the fitness 
cost and the fungicide concentration 



Prb = k k /(k k + 1), 



(A.21) 



C b = C 50 / tz P \ • (A.22) 

«k(l - Pr) ~ Pr 

In this case the sensitive strain dominates at C < if p T < p r b or at any positive values 
of C if p r > p r b [white area in Fig.[2^b)]. 

When only the low risk fungicide (fungicide B) is applied, we set r-Q = 1 in the 
inequality (1A.8I) and obtain 

<Pr/k k + - % . (A.23) 



C + C50 C + C50 



This inequality holds and the sensitive strain dominates for all positive values of p r and 
C at which R 0s > 1 

Consider the case when the two fungicides A and B are applied together at an arbitrary 
mixing ratio re, assuming no pharmacological interaction (u = 0) and perfect resistance 
of the resistant strain to the fungicide A (a = 0). In this case, 7„ = 1 and 7 r = 
r B . Substituting these values in Eq. (IA. 13[) . Eq. flA. 14[) and Eq. flA. 15[) gives the same 



expressions as in Eq. flA.18j) . Eq. flA. 19[) and Eq. flA.20p . but with a substituted by r B . 



A. 3. Expressions for the treatment benefit 

The treatment benefit is defined as the ratio between the amount of susceptible hosts 
H(t) when both the disease and treatment are present and the amount of healthy hosts 
at no disease B(t) = H(t)/H nd (t) (see Sec. l3.3p . We consider here its equilibrium value 
over long time scales given by 

H* 

B{t -> oo) = B* = — , (A.24) 
K 

where H* is the equilibrium amount of healthy hosts and K is the host carrying capacity, 
and assume full resistance (a = 0). 
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When only the fungicide B is applied at a concentration C [Fig.H^a)], the basic repro- 
ductive number of the sensitive pathogen strain always exceeds the one for the resistant 
strain R 0s > i? 0r . Therefore, the resistant mutants are eliminated in the long run and 
the amount of the healthy host tissue is equal to H* = fi [1 + e(C)], where e(C) is given 
by Eq. fl4]). Then, according to Eq. flA.24jl . the treatment benefit is 

B*(C)= p[l + e(C)}/(bK). (A.25) 

It grows with the fungicide concentration and saturates, since the function e(C) satu- 
rates. 

Application of the fungicide A alone at a concentration C may favor either resistant or 
sensitive pathogen strain depending on the fitness cost of resistance p r and the fungicide 
concentration C [see Fig.^b)]. The treatment benefit in this case is 

+ e(C)] /(bK), for (p r < k k /(l + k k ) and C < C h ) 
B*(C, p T ) = I or (p r > k k / {I + k k ) and VC), (A.26) 

l/V [b(l - p r )K] , for p r < k k /(l + k k ) and C > C b , 



where the Cb is given by Eq. flA.22j) . 



Now, consider application of both fungicides in a mixture at equal concentrations 
(tb = 1/2); assuming no interaction between fungicides (u = 0). In this case, again 
either resistant or sensitive pathogen strain will dominate the population depending on 
the fitness cost p r and the total fungicide concentration C [see Fig.|5(c)]. The treatment 
benefit now has the following expression 

U [1 + e{C)\ /(bK), for (p r < p rb and (C < C hl or C > C h2 )) 

B*(C,p T ) = < or (p r > p rb and VC), 

[fi[l + e(C/2)} I [6(1 - p r )K] , for p T < p rh and C bl < C < C h2l 

(A.27) 



where p r b, C b i and Cb2 and are given by Eq. (1A.13I) and Eq. (1A.14j) at 7 S = 1, 7 r = 1/2. 



The boundary value of the resistance cost has then a simpler expression p r b = k^(2kk 
2y/2(l + kk) + 3)/(l + 2k k f. The treatment benefit B*(C,p r ) is shown function of 
the fungicide concentration C and the fitness cost of resistance p r in Fig. H] for the three 
cases discussed above, according to Eqs. (]A.25ft - (]A.27|) . 



A. 4. Selection for resistance is delayed by applying a fungicide 
mixture 

If the fungicide resistance is not associated with a fitness cost, then the resistant strain 
is favored by selection and eventually dominates the population whenever the high risk 
fungicide is applied alone or in a mixture with the low risk fungicide [Fig.[2£b,c)]. How- 
ever, for a given value of the total fungicide concentration C, the selection for resistance 
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slows down when applying the fungicide mixture as compared to the treatment with the 
high risk fungicide alone [as seen from time- dependent numerical solutions of the model 
Eqs. ([I])-©] in agreement with the findings of [3]. 

In order to understand this result we consider the dynamics of the frequency of the 
resistant pathogen strain p{t) = I T /(I r + J B ). The rate of its change is obtained from 
Eqs. ©-© [9] 

f t =sp(l-p), (A.28) 

where 

s = fi [e s (C, r B ) - e T (C, r B )] - p T bH(t) (A.29) 

is the selection coefficient [a similar expression was found in [50]]. Here e s (C, r B ) and 
£r(C, r B) ar e given by Eqs. (j5J, (EJ) and r B is the proportion of the fungicide B in the 
mixture. 

Assuming a zero fitness cost (p r = 0), no pharmacological interaction (u = 0) and full 
resistance (a = 0), the selection coefficient can be written as 

s = ii [e{C) - e(r B C)] . (A.30) 



In this case, the selection coefficient does not depend on time and the Eq. (1A.28|) has a 
simple analytical solution: 

Pit) = 1 P °t MSt \ TT ' (A-31) 

where po = p(t = 0). At s > 0, the function p(t) grows monotonously and tends to 
one at large times. The rate, at which it grows is determined by the magnitude of the 
selection coefficient s. 

One can see from Eq. (1A.30I) that when the high risk fungicide is applied alone (r B = 0), 
the selection coefficient is larger than when it is mixed with a low risk fungicide (0 < r B < 
1) at the same total fungicide concentration C. Hence, s(r B = 0, C) > s(r B > 0, C). 
This is because the function e(r B C) has positive values for any r B > 0. Thus, the 
selection for the resistant strain (against the sensitive strain) is delayed when a mixture 
of high risk and low risk fungicides is applied compared to treatment with the high risk 
fungicide alone. 



A. 5. Generalization of the model 

So far we assumed that resistance cost manifests in the transmission rate b of the resistant 
pathogen strain and the fungicide increases the pathogen death rate /i. We performed 
the same analysis for the three remaining cases possible in the model: When (i) both 
resistance cost and the fungicide affect the pathogen death rate according to \i — > + 
p T +e r (C, r B )) for the resistant strain and \x /i.(l+e s (C, r B )) for the sensitive strain; (ii) 
both resistance cost and the fungicide affect the infection rate b — > b(l — p T — e r (C, r B )) for 
the resistant strain and b — > b(l — e s (C, r B )) for the sensitive strain; (iii) resistance cost 
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affects the death rate of the resistant pathogen strain fi — > + p r ), while the fungicide 
affects the infection rate of both resistant and sensitive strains b —> b(l — e T)B (C, vb)). We 
have found that although the mathematical expressions for the results have a different 
form in these cases, all the conclusions remain the same, so that the results do not 
depend on whether the fungicide and the resistance cost manifest in the infection rate 
b or in the pathogen death rate \x. Moreover, we have done the same analysis using 
a fungicide dose- response function different from Eq. (j3J), namely using an exponential 
function and again found that all the conclusions remain the same. 

This generalization applies to the results at equilibrium, but the time-dependent so- 
lutions of Eqs. ([I])-© may behave differently depending on how the fungicide and the 
fitness cost affect the pathogen life cycle and the form of the fungicide dose-response 
function. This is an interesting topic for further investigations, but lies beyond the scope 
of this study. 

A. 6. Fungicide mixture versus alternation 

It was previously discussed [15] that in the presence of a fitness cost the alternation of 
fungicides can be effective, but we have shown here that fungicide mixtures can also 
be effective in this case. When using an alternation strategy, the period of selection 
during which the resistant strain is favored in the presence of the high risk fungicide is 
followed by a period during which selection favors the sensitive strain in the absence of 
this fungicide. The latter period is typically much longer because the selection pressure 
induced by the high risk fungicide is much larger than that induced by the fitness cost 
of resistance. Hence, one needs to wait for quite a long time before the resistant strain 
disappears and the high risk fungicide can be used again. Moreover, there are times 
during which the frequency of the resistant strain becomes large (at the end of the period 
of the application of the high risk fungicide), which increases the risk that resistance 
will spread to other regions. Both of these disadvantages are avoided by using a mixture 
where the proportion of the low risk fungicide is above a critical value determined here 
(Fig.[3]). In this case there is no need to delay the application of the high risk fungicide 
and the frequency of the resistant strain does not rise above the mutation- or migration- 
selection equilibrium because the mixture does not induce selection for resistance. 

A. 7. The risk of double resistance 

Although we do not consider the possibility of double resistance in our model, by ap- 
plying an optimal proportion of fungicides in the mixture as suggested here, one would 
prevent selection for resistance to the high risk fungicide. Consequently, the risk of 
development of double resistance would be reduced. For both sexually and asexually 
reproducing pathogens, there are three pathways for generating double resistance: (i) 
A-resistant mutants are produced first and then a proportion of them acquires also 
B-resistance by spontaneous mutation (ii) B-resistant mutants are generated first and 
subsequently acquire A-resistance and (iii) double resistance is generated directly from 
the wild-type. In this case, by preventing selection for A-resistance, one removes only 
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the pathway (i) to double resistance. If a pathogen is able to reproduce sexually, then 
a much more likely scenario for the double resistance to emerge is through recombi- 
nation. For the recombination to occur, both singly resistant strains (A-resistant and 
B-resistant) would need to be present in the population at significant frequencies. Hence, 
preventing selection for A-resistance would diminish the probability of the emergence of 
double resistance by recombination. Thus, our findings would also help to significantly 
reduce the risk of development of double resistance, especially in sexually reproducing 
pathogens. 
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